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Abstract 

We have studied the low-energy excitation spectrum of a dimerized and frus- 
trated antiferromagnetic Heisenberg chain. We use an analytic approach, 
based on a description of the excitations as triplets above a strong-coupling 
singlet ground state. The quasiparticle spectrum is calculated by treating the 
excitations as a dilute Bose gas with infinite on-site repulsion. Additional 
singlet (S=0) and triplet (S=l) modes are found as two-particle bound states 
of the elementary triplets. We have also calculated the contributions of the el- 
ementary and collective excitations into the spin structure factor. Our results 
are in excellent agreement with exact diagonalizations and dimer series expan- 
sions data as long as the dimerization parameter 6 is not too small {5 > 0.1), 
i.e. while the elementary triplets can be treated as localized objects. 

PACS: 75.10.Jm, 75.30.Ds, 75.40.Gb 
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I. INTRODUCTION 



The properties of a variety of recently discovered quasi one-dimensional materials can 
be described by the Alternating Heisenberg Chain (AHC) model. The Hamiltonian of the 
model reads: 

H = JY^[{1 + 5(-l)0S,.S,,+i + «S,.S,+2], (1) 

i 

where i denotes the sites of a chain with length and Sj are S = 1/2 spin operators. 
The parameter a is the next-nearest neighbor coupling, leading to frustration, and 5 is the 
dimerization. We assume J > 0. 

The AHC is realized in nature in many materials that have two important but structurally 
inequivalent superexchange paths that are spatially linked, so that an alternating spin-spin 
interaction results. Representative examples of materials of this type are {yO)2P20-j |]l|] and 
various aromatic free-radical compounds 0. Alternating chains may also result 
of spontaneous dimerization, due to dynamical spin-phonon coupling. Examples include the 
recently discovered inorganic spin-Peierls compounds CuGeO^ [0] and a' — NaV20^ 0]. 

The AHC model (|l]) is a straightforward generalization of the uniform Heisenberg chain, 
which is the most widely studied quantum spin system. The uniform 5 = 1/2 chain has 
a gapless excitation spectrum with a known dispersion relation and a rather complicated 
ground state which is characterized by strong quantum fluctuations, making it highly unsta- 
ble to perturbations. The AHC generalizes the uniform chain by alternating the spin-spin 
interaction between two values, J(l + 5) and J(l — 5), and including next-nearest neighbor 
coupling (frustration) aJ. Since the Hamiltonian (|l|) is rotationally invariant with respect 
to spin, the total spin is a good quantum number, and the ground state is a spin singlet. The 
translational symmetry however is broken by the dimerization, and the resulting system has 
a gap to the first excited state with S = 1. The model has a rather complicated spectrum 
of states at higher energies, including multimagnon continua and bound states. 

Recent neutron scattering studies of {yO)2P207 carried out by Garrett et al. 
show that this material can be described well by the alternating spin chain model. In 
addition, these measurements provide evidence for the existence of a two-magnon triplet 
[S = 1) bound state. Similar studies of the spin-Peierls material CuGeO^ by Ain et al. [0, 
in combination with Raman spectroscopy p|, suggest that a singlet 5* = bound state is 
present in this material. These experiments have strongly motivated the theoretical studies 
of the low-energy excitations of the AHC model, since the existence of magnetic bound 
states is one of its characteristic features. More generally, bound states seem to appear 
in all low-dimensional dimerized quantum antiferromagnets which have a gapped triplet 
spectrum above a singlet ground state [§. 

Besides the relevance to real compounds, the model (p is interesting from theoretical 
point of view since it contains two independent mechanisms for spin gap formation. Subject 
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of a special interest is the nature of excitations. At 5 = the model exhibits two phases, 
separated by a critical value of frustration ac = 0.2411, known from numerous studies 
p!0| . For a < Oc the ground state is similar to that of the uniform Heisenberg chain and 
frustration is irrelevant in this parameter regime. The elementary excitations are massless 
unconfined spinous with spin S=l/2 [|lT|. At a = ac there is a transition into a phase with 
spontaneously dimerized two-fold degenerate ground state. The spectrum acquires a gap 



and the elementary excitations are massive spinous |1T2|-[T^. On the other hand Haldane |T2 



has shown that for any 6^0 the spinous become confined into triplet S = 1 excitations 
with a gap in the spectrum. Interactions between the triplets can lead to additional massive 
singlet, S = 0, and triplet, S = 1, excitations below the continuum |]15|,|13| and even a 
sequence of further massive excitations [|1^]. As 5 — the size of the elementary triplet (or, 
equivalently, the spinon radius of confinement), Tq, is known to scale as ||T7| , p^ : Tq ~ 5^^/^ 
for a < ac and tq ~ S^^^^ for a > ac- Therefore the triplets become nonlocal objects for 
small dimerization, ultimately giving way to completely unconfined spinous in the strict 
hmit 6 = 0. 

The excitation spectrum of the AHC has been extensively studied by a variety of nu- 
merical methods. Most recent work includes: Exact Diagonalizations (ED) carried out by 
Bouzerar, Kampf, and Japaridze [jl9[ and Barnes, Riera and Tennant [BO], dimer series ex- 



pansions by Singh and Weihong [^, and density matrix renormalization group (DMRG) 
studies by Sorensen et al. p2|. However there are only few analytical works studying the 
AHC model. Uhrig and Schulz |]15| have used the continuum limit field theory approach and 
predicted the possibility of singlet (S=0) and triplet (S=l) two-magnon bound states below 



the two-particle continuum in the dimerized phase (5 > 0). Affleck and Uhrig et al. |17 



developed a description, based on the soliton picture in which the elementary excitations are 
treated as confined 5 = 1/2 solitons (spinous). It was demonstrated that = 1,0 soliton 
anti-soliton bound states are formed and their number increases as 5 0. Bouzerar and Sil 
p3| have studied the excitation spectrum of the dimerized phase using the Bond Operator 



Technique (BOT) introduced by Chubukov and Sachdev and Bhatt p5|, in combination 
with the Brueckner approach developed by Kotov, Sushkov, Weihong, and Oitmaa |^|. The 
results of Ref. p3[ are mainly inconsistent with the ED data. They are also inconsistent 
with our results in spite of the fact that we use the same approach. The reason for this 
discrepancy, as will become clear from our presentation below, is that the treatment of Ref. 
p3| is quite incomplete and misses several important technical aspects of the problem. 

The aim of the present work is to investigate theoretically the excited states of the 
dimerized phase {5 > 0) of the AHC by using the analytic Brueckner approach. Previous 
applications of the Brueckner approach include the two-layer Heisenberg model p6| , p7| , the 
quantum spin-ladder as well as the spin ladder with frustration ||2^ . In all cases excellent 
agreement was achieved between the theoretical and numerical spectra. The Brueckner 
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approach is based on the description of the excitations as triplets above a strong couphng 
singlet ground state. In essence it represents an effective way of taking into account the 
hard-core constraint, which has to be imposed on the triplets, by treating the excitations as 
a dilute Bose gas with infinite on-site repulsion. The method is valid while the on-site density 
of excitations (triplets) is small and under the assumption that the elementary triplets 
are well localized objects. In spite of the fact that the on-site density in the AHC is small 
even close to the point 6 = [rii ^ 0.08 for 6 = 0.01, a = 0), we find that our approach fails 
at (5 ~ 0.1 because the triplets become nonlocal objects. In the region 6 > 0.1 the present 
method provides an excellent quantitative description of the excitation spectrum. 

The rest of the paper is organized as follows. In Sec. II we find the system of self-consistent 
equations describing the AHC with frustration and calculate the one-particle (elementary 
triplet) spectrum for a wide range of parameters. Sec. Ill describes the additional singlet and 
triplet modes which appear below the two-magnon continuum. These excitations are bound 
states of two elementary triplets. In Sec. IV we calculate the contributions of the elementary 
and the additional triplet into the neutron scattering structure factor. Sec.V contains our 
conclusions. 

II. FORMULATION OF THE PROBLEM AND QUASIPARTICLE SPECTRUM. 

In order to analyze the excitation spectrum of (|1|) it is convenient to adopt the strong- 
coupling viewpoint. In the limit 6 = l,a = the ground state consists of non-overlapping 
spin singlets \GS >= |1,0 > |2,0 > |3,0 > where \i,0 >= -^[\ t>2i | i>2i+i -| i>i 
I T>2i+i]- Here, and in all formulas below which involve triplet operators, the "generalized" 
site index i represents the bonds with spin exchange J(l + 6). Each singlet can be excited 
into a triplet state and therefore it is natural to introduce a creation operator for this 
excitation: \i, a >= tl^i\i, >, a = x,y, z. The representation of the spin operators in terms 
of the triplets was introduced in and reads: 



Sa,2i — 2^^"* ^ai '^(^af^-ytpit^yi) , (2) 



2^ -ui -trt -taP'fl^fSi 

where eaf^^ is the fully antisymmetric tensor. After application of this transformation to (p, 
or, equivalently, after calculating the matrix elements of the "hopping" terms, we find the 
effective Hamiltonian 

H = H2 + H; + H^ + Hu, (3) 

N/2 ^ N/2 

H2 = J±Y. titai + 2 E(4Wi + 44+1 + h.c), (4) 

i=l,a ai 
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j=l,a/37 



i=l,/37 

Hu = i7 44^/3^^"^' ^ 00, (7) 

^ i=l,a/3 

where we have adopted the notation: 

J^ = J{l + 5),\ = -^^(1 - 5 - 2a), /i = i J(l - (5 + 2a), z/ = J(l - 5). (8) 



The present notation for the coefficients is chosen following Ref. [28|. For the problem 



studied in Ref. the Hamiltonian represents a spin ladder with frustration and the term 
cancels out due to the symmetry of the ladder with respect to interchange of its two 
legs. The coupling J(l + 5) corresponds to the rung exchange Jj_ in the spin ladder model, 
and gives the on-site energy of the triplets (ffist term in Eq.(^). The rest of the terms in 
the Hamiltonian represent an effective hopping of the triplets, spontaneous creation of two 
triplets on nearby sites (second term in Eq.(|)), as well as interactions between the triplets 
(Eq.(||) and Eq.(|^)). In addition, an infinite on-site repulsion Hu is introduced to take into 
account the hard-core constraint t^jt^j = 0, which has to be satisfied on every site. This 
condition (only one triplet can be excited on a site) follows from quantization of spin and 
guarantees the one to one correspondence between the spin model Eq.(|]) and its effective 
triplet description Eq.(|^). 

Let us ffist consider the Hamiltonian (|^) without the i^a term. A detailed analysis 
of the spectrum was presented in Refs. p6|,p|,p8| an we will only summarize the results 
here. The interaction Hjj gives the dominant contribution to the renormalization of the 
triplet spectrum. It was demonstrated in [26] that in the Brueckner approximation this 



renormalization is given by the normal self-energy operator 

^n{k.^) = 4E ^'r^'^(A; + q,u- u,), (9) 

where the scattering amplitude is 

/ r^Brr^Br 2 2 \ 

The renormalized spin-wave spectrum Uk, the quasiparticle residue Z^. and the Bogoliubov 
parameters Uk, Vk can be found from the solution of the coupled Dyson equations for the nor- 
mal Gn{k,t) = -i< T{tka{t)tl^{0)) > and anomalous Ga{k,t) = -i < T(tl;.„(t)4„(0)) > 
Green's functions (see Ref. IEgI): 



5 



^k = Z^^^Al-Bl (11) 

ifc = + Acos2/c + T.^"" {k, 0) + 2/icos2A; J2 Zg'v^^cos2q, (12) 
-Bfc = Acos2A; — 2/icos2A; ^ Z^^UqVqCos2q, 

^f=(l-^ ) > (13) 

2 2 

These equations also take into account the quartic interaction, (||), in the lowest order, 
one-loop approximation. For the sake of simplicity, in all of the above formulas, and every- 
where from now on, we omit the number of lattice sites. All momenta take values inside 
the Brillouin zone of the original (non-dimerized) lattice (— vr/2 < k < tt/2). The discrete 
lattice sums can be converted into integrals by J2k 

The basic approximation made in the derivation of r^'^{q,uj), as discussed in Ref. p6 




is the neglect of all anomalous scattering vertices, which are present in the theory due to 
the existence of anomalous Green's functions (i.e. quantum fluctuations). All anomalous 
contributions however are suppressed by the small parameter of the Brueckner approach: 
the density of triplet excitations rii = J2a < tlatia >= ^^Y^q Z^'^Vg. As follows from Eq.(|D 
the Brueckner self-energy is first order in Ui, while it is easy to prove that all anomalous 
contributions lead to higher powers of the density This observation justifies the Brueck- 



ner approximation (i.e. keeping only Eq.(P)), since the triplet density is small throughout 
the disordered phase, even close to 5 = 0. We find that rii ^ 0.0066 for 6 = 0.6, a = 0, 
increasing at decreasing dimerizetion rii ^ 0.066 for 5 = 0.1, a = 0. 

Now consider the effect of the cubic interaction H3, Eq.(^). In order to incorporate the 
effect of self-consistently, the normal and anomalous self-energies have to be found to a 
given order and subsequently inserted into the system of coupled Dyson's equations for the 
anomalous and normal Green's functions. The poles of these Green's functions determine 
the renormalized one-particle spectrum [Q. This procedure however is rather involved and 
in the present work we will use only the lowest order perturbation theory result (i.e. we will 
keep one diagram only). This is justified, since we find that the lowest order contribution is 
relatively small and therefore does not require a fully self-consistent treatment. The cubic 
interaction Eq.(^, after Fourier: taj = Hk'^ake^^^^ ■, and Bogoliubov: tka = UkCika + Vka}_ka 
transformations can be rewritten as 

Hi = - Yl ^"/37 (ri(l, 2, 3)[al^„al2^afe3^ + h.c] + (14) 

ki,k2,k3=ki+k2 
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+ lr2(l,2,3)[a[ 



(15) 



We have defined 



ri(l, 2, 3) = (Mfc3Z^(l, 2, 3) - Vk,C{l, 2, 3)) 



7Br 



yBr 
^k2 



yBr 
^k3 ' 



C(l, 2, 3) = ro(l, 2)vk^Vk2 + ro(-3, 2)uk2Vk^ + ro(l, -3)nA:iffc2, 

where ro(l,2) = — z/sin(fci — A;2)cos(A;i+/c2) andD(l,2, 3) = C(1,2,3) with replacing all Uk by 
t>k and vice versa {u ^ v}, and r2(l, 2, 3) = — ri(l, 2, 3) with replacing {uk^ ^ ^k^}. Only a 
normal Green's function Gn = —i < T{aka(t)(4cai^)) > exists for the physical operators a^a, 
and therefore the self-energy induced by is given, to lowest order, by the two diagrams 
shown in Fig.l: 



'Ti{q,k-q,k) Tl{q,k-q,ky 



UJ — UJn 



^k- 



(16) 



The correction 6uk to the one-particle spectrum is: 6u!k 
also gives a correction to the quasiparticle residue: zj?^ ; 
In order to find the spectrum, equations 



T,3{k,ujk). The self-energy (|T6D 



1 + 



du) 



|0],^|12],|13D have to be solved self-consistently 
for S^'"(fc, 0) and Z^^', which corresponds to an infinite re-summation of diagrams (solution 
of Dyson's equation), and leads to the renormalized spectrum ojk- Then the correction due 
to if 3 has to be added resulting in the spectrum VLk = oJk + SuJk- In this approximation the 
quasiparticle residue is Zk = Z^^zj^K In the limit z/, A ^ J_l all expressions can be easily 
evaluated analytically to leading order, and for example the corrections, due to H-^ are: 
6uJk = — ^j-cos^/c, and zj^^ = 1 — ^4-cos^fc. Finally, combining the three-particle corrections 



with the results of Ref. for the one-particle spectrum without H^, we obtain in the limit 
A, 1/ < Jr. 



Qk 

j7 



A , 3 A2 
1 + — cos2fc + -— 5" 



1^' 2, 
TrCOSAk ^cos k. 

4 Jf 2 Jf 



(17) 



Higher orders also can be evaluated, similarly to Ref. but we shall not present the 
lengthy expressions here. 

Next we present results obtained by numerical self-consistent solution of the Dyson's 
equations. Plots of the triplet gap A = Qk=o/J as a function of a at 5 = 0.4,0.2 are 



presented in Fig. 2. Our results agree quite well with the ED data |T9[. We also note that 



our results disagree with those presented in Ref. The main reason for this is that the 
three-leg vertices H^, Eq.(|D were not taken into account in Ref. [^. In figures 4a,b,c,d 
we also present spectra of the elementary triplet Qk for different parameter values. The 
agreement with the numerical data (when available) is very good. 

The approach presented here relies heavily on the effective description of the spin problem 
in terms of localized triplets. As mentioned in the Introduction, in the limit 6 —>■ 0, for any 
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a, the spinons, confined to form tlie triplets, become unconfined. Our calculation clearly 
breaks down at this point. We would like to emphasize however, that due to the power law 
dependence of the radius of confinement tq on 5, namely tq ~ for a < = 0.2411 

and tq ~ 5~^^^ for a > ac [pT| , pTl , the triplets are well localized for S > 0.1. For S ~ 0.1 the 



triplet size is about 2-3 lattice spacings. Therefore the strong-coupling approach in fact is 
expected to describe quantitatively well a very wide range of parameter space. 

Let us also mention that the effective Hamiltonian Eq.(^ can be used even for 6 < 0.1, 
when ro increases strongly. However technically the problem becomes much more difficult 
and truly non-perturbative in nature. One way to deal with it is to develop a perturbation 
theory to high order, e.g. via the dimer series expansion ||21|]. Physically a finite ro leads 



to the creation of a "string" (with length tq) of excited triplets. This string leads to the 
formation of low-energy many-particle bound states which mix strongly with the one-particle 
excitations, causing ultimately (for 5 — 0) the vanishing of the one-particle spectral weight. 
These effects can be taken into account by diagonalizing the Hamiltonian in Fock space 
(spanned by the one-magnon, two-magnon, etc. states). A detailed analysis however is 



beyond the scope of the present work |31 



III. SPECTRUM OF COLLECTIVE EXCITATIONS: TWO-MAGNON BOUND 

STATES. 

The quartic term Eq.(^ in the Hamiltonian Eq.(|^) leads to attraction between two 
elementary triplets. It has been already demonstrated in the case of a spin ladder PJ2E[| 
that this attraction is strong enough to form additional singlet (S=0) and triplet (S=l) 
bound states below the two-particle continuum E^{Q) = m.mg{flQ_^ + nQ^^). Here we will 
calculate the energies of these additional singlet and triplet states in the AHC model. Such 
bound states have already been observed in the ED calculations P^,|2D[] of the AHC. 



Consider the scattering of the two triplets: qia + g2/3 qs'J + q^S and introduce the 
total, Q, and relative, q, momenta of the pair. The two-particle singlet and triplet wave 
functions are: 

l^a(Q) >= ^ E ^^iQ,QKp,a^Q^yQ_jO >, (19) 

where ip^''^{Q, q) are determined from the Schrodinger equation, H\ip^''^ {Q) >= 
E^''^(Q)|V^^'^(Q) > with bound state energy E^''^{Q). From this equation one can read- 
ily derive the integral Bethe-Salpeter equation satisfied by the bound state wave-functions: 

[E'^'^iQ) - nc^. - n<iU''^{Q,q) = J2M''^{Q,q,p)^''^{Q,p). (20) 
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Let us at first neglect the cubic interaction i^s in the Hamiltonian (^. Then the scattering 
amphtudes M^''^{Q, q,p) in the singlet and triplet channels are (see Ref. |2^): M^{Q, q,p) = 
(f/ — 2/xcos2gcos2p), and M^{Q, q,p) = — yusin2gsin2p, respectively. Eq.(pOD should be solved 
with the substitution 

M^'^(Q,g,p) ^ y^y^y^y^n^_^«^^^«^_^«^_,^M^'^(g,g^ (21) 

in order to take into account that the triplet excitation al^ differs from the bare one 4^,, 
due to the Bogoliubov transformation and the quasiparticle residue. We stress that the 
solution of the integral equation (EO) in the singlet channel has to satisfy the condition 



Y.qi^^ {Qil) = 0, due to the infinite on-site repulsion. A Lagrange multiplier has to be 
introduced to enforce this condition. The normalization constants of the wave functions 
(|I^,|T^) are chosen to satisfy the conditions: J2q |V'(Q)P = J2q \4'{Q,(l)\'^ = 1- The solution 



of Eq.(pOD for the triplet and singlet bound states (0,|1^) has already been found in [^] for 
the Hamiltonian (§) without the if 3 term. In the leading order in X/ J±, fi/ J± the energies 
and the wave functions of the bound states are: 

E'q = 2J^ - /.(I + C'q) + i^^os2Q, (22) 



^^(Q,g) = ^2(l-C2)- 



cos2g + Cr 



2CQCos2g + 



2 ' 



= 2J^ - ^(1 + 4CJ) + If + -Vcos2Q, (23) 



sin2g 



2 +2CQCOs2g + 2C2 ' 



where we have defined Cq = ^cosQ. 

Next, we analyze how affects the two-particle spectra. The contribution of the three- 
particle scattering into binding should change the form of Eq. (^) since it leads to retardation 
and thus to a non-trivial frequency dependence in the Bethe-Salpeter equation. In the 
present work we treat as a perturbation and find the second order corrections, 6Eq , to 

S T 

the two-particle bound state energies Eq . The correction in the singlet channel is given as 
a sum of the diagrams, shown in Fig.3a,b,c, convoluted with the singlet wave function 

5Ef{Q) =Y^5M^{Q,q,p)^''{Q,q)^'{Q,p) + 5E^{Q) + 5E^{Q), (24) 
with the vertex 

SM^(Q, ,.p) = _ 4r.(|-,,,-p.f-p)r.(f + p.,-p.a + ,) 
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and 



5E!{Q) = 2Y.rsiQ,p) J^^^'^ • (27) 

The correction to the triplet energy E^{Q) is given by the diagrams a,b,c,d,e in Fig.3, 
convoluted with the triplet wave function 



6El{Q) = J2{6Mj + 6Mj + 6MJ)^''{Q, q)V{Q,p) + SE^iQ) + SEJ{Q) (28) 

p,q 

with the vertices: 5Mj(Q,g,p) = i(5Mf (Q, g,p){E^(g) E^{Q)}, 
5MjiQ,q,p) 



ri(f +P,f -p,Q)ri(f + g,f -g,Q) 



SMj{Q,q,p) 



r2(§ + g,§-g,Q)r2(f + p,f -P,Q) 



r+9 ^f-p ^f+p 
and the terms: SE^Q) = 6E^^^{Q){^^, E^ E^}. 

In the leading order in t^, A ^ J± the corrections (0,^) can be easily found analytically: 

SE^iQ) = -^cos'Q, SEfiQ) = -^(1 - ^cos^Q). (29) 

We observe that the contributions (|2^J2^) have an appreciable effect on the two-particle 
energies and therefore the cubic interaction (1^) can not be neglected. 



It is worth noting that Eqs. (PO j^ psD are not quite correct because they do not take fully 
into account the hard-core constraint tj^tj^ = 0. Let us consider the limit X,!/ <^ fi <ti J±. 
A direct expansion in powers of 1/J± applied to the Hamiltonian (|^) in coordinate space, 
gives the following energies for the singlet and triplet bound states: 

Q\2 ^2 

ELi = - - ^(1 + cos2g) + ^ - ^(1 + COS2Q), (30) 

,, \2 q\2 2 

ELi = 2Jx - f - -(1 + COS2Q) + ^ + --COS2Q. (31) 
Z jj, oJ± oJ± 

The hard-core constraint is completely taken into account in these formulas by not allowing 
hopping processes leading to two triplets occupying the same site. On the other hand the 



momentum space diagrammatic calculation Eqs.( P2| , P3D and the three-particle correction 
(p9|) result in 

A^ A^ z/^ 

= 2Jx - - —(1 + COS2Q) + ^(6 - cos2g) - ^cos^Q, (32) 

= 2J^ - I - -(1 + COS2Q) + ^(6 + cos2g) " ^(1 - 2^08^^). (33) 
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The physical reason for the difference between the momentum space calculation and the exact 
real space result is in the additional blocking of virtual excitations which is not included 
in Eqs. ( pO| , p^p8D . Let us consider one example of such additional blocking. When we 
have a state with one triplet, the quasiparticle blocks virtual excitations due to quantum 
fluctuations, tJa^f+ia' t"^*^ links. This is the physical origin of the third term in the 
one-particle dispersion (p!?!). When we have two excited triplets and they are separated by 
more than one lattice spacing they block four links (this corresponds to the third terms 
in (p2| , p3D ). However when the two triplets are on nearest-neighbor sites they block only 



three links, which is energetically more favorable and leads to the increase of the binding 
energy. This effect represents a non-potential contribution and is not taken into account 
by the integral equation (^). The additional contribution of the quantum fluctuations into 
binding is analogous to the calculation of the Lamb shift in atomic physics. It is clear from 
the above discussion that binding due to blocking of quantum fluctuations exists only when 
the triplets are on nearest-neighbor sites. 

We have also identified several additional processes which have to be described more 
accurately by inserting the vertex (|10[), responsible for the infinite one-site repulsion, into 
the diagrams presented in Fig. 3. These diagrams should account correctly for the hard-core 
repulsion between the particles in the initial and intermediate states. However it is very 
difficult to identify and calculate in momentum space all diagrams of this type. We will follow 
a simpler, effective way to take into account these additional contributions. Comparing the 
coordinate space results (|30|j31| ) with the momentum space calculations (|3^ , |33|) one can see 
that the latter can be corrected by adding the following expressions to Eqs.p^j28|): 



6E^{Q) 
6E^{Q) 















8J± 





3 
2 



cos2(5 
2 

-cos2( 
3 



^v^cos2A;^^(Q,A;) 

^y2sin2A;V^^(Q,A;) 



A2 



(34) 
(35) 



In the above equations the expressions in the absolute value signs give the probability am- 
plitudes for two quasiparticles to be on nearest-neighbor sites in the appropriate channel. 
Finally, to obtain the energies of the triplet and singlet bound states, the integral equation 
(13) should be solved first, and then the corrections SE^'^{Q) (||,|28|) and SE^''^{Q) (HH) 



should be added. Results obtained by a self-consistent numerical solution for the spectrum 
of the singlet and triplet bound states are presented in Fig.4 for different values of the frus- 
tration, a, and dimerization, 6. We find that there are no bound states in the vicinity of 
q = for zero frustration (a = 0). Singlet and triplet bound states are always present in 
the vicinity of g = 7r/2, however the singlet typically exists within a much wider momentum 
range. At g = the singlet splits off from the continuum at any non-zero a, while the triplet 
exists only above a certain value of frustration. One can see also that, quite generally, the 
singlet is below the triplet, i.e. the binding in the singlet channel is stronger. 
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We have also calculated the binding energies, e^''^{Q) = E'~^{Q) — E^''^{Q), of the singlet 
and triplet bound states as a function of 5 at Q = 7r/2 for a = 0. The results are presented 
in Fig. 5. There is excellent agreement with the ED data [p!9| , p0| in the singlet channel even 
for 6 = 0.1, whereas the agreement for the triplet is not so good for 6 < 0.5. We attribute 
this disagreement to the fact that the triplet binding energy is relatively small in comparison 
with the singlet one for the whole range of parameters {a, S). Therefore taking into account 
the three-particle scattering in simple perturbation theory is not as good approximation for 
the triplet as it is for the singlet. We have performed a more accurate calculation for the 
cases when the binding energy is very small. An improvement can be achieved by adding the 
corrections SM^''^(Q, q,p) to the scattering amplitude M^''^{Q, q,p) into the Bethe-Salpeter 
integral equation (pO]) and further solving the equation to find the binding energies. Such 
"self-consistent" calculations are presented by the dashed lines in Fig. 5 and one can see that 
the agreement with the ED data indeed becomes better for small energies. 

There are several important points which have been overlooked in Ref. ||2^, ultimately 
leading to incorrect results for the energies of the bound states. The contribution of the 
constraint has not been taken into account in the integral equation (^) resulting in 
non-zero binding energy of the singlet at g = even for large S and a = (notice that our 
calculation always gives zero binding at this point). The contribution of the three-particle 
scattering (^) into binding has not been taken into account in Ref. ||23| and the effect given 
by the corrections ( p4| , |35D (Lamb shift) has not been included. In addition, the error in the 



calculation of the one-particle spectrum (see the previous section) has propagated into the 
two-particle energies as well. 

Finally we would like to compare our results for the bound state energies with the lowest 
order perturbation theory results, presented in Ref. The binding energy e^''^{Q) = 

E'^{Q)-E''^{Q)atQ 
(1-6) 



n/2, a 



SfTT 



e 



14(1-5)2 



17 (i-sy 

32 l+<5 



resulting from (|30| , |3TD and ([T^) in the singlet channel. 



is different from the result of Ref. 



01: e'{l)/J=h 



g2 • In the triplet channel our result is 6e^{^)/J — ^ 
with the formula presented in Ref. [pO . 



^^1 — 6) — II , which coincides 



IV. STRUCTURE FACTOR. 



The one-particle triplet and the triplet bound state can be observed in neutron scattering 
experiments. The inelastic neutron scattering cross-section is directly proportional to the 
dynamic structure factor: 

Siq, u) = hm J e^-* < T(S(g, t).S(-g, 0)) > . (36) 

Here we calculate the spectral weights of the elementary and the additional triplets. After 
the bond operator (^ and Bogoliubov transformations, the on-site spin operator Sj can be 
written in momentum space as: 
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S^{-q,0) = -te"'/' 



(37) 



where we have defined the vertices Ti(g) = {ug + fg)sin| and T2{k,q) = Uk+q/2Vk-q/2COSj. 
Then, after averaging, the structure factor can be rewritten in the form S{q, uj) = Si{q)6{uj — 
+ S2{q)6{u; - Ej) where 



Si{q) 



'^T. (,) + 2 V T.ik,q)T,ik + l^-k + l,q) 



(38) 



is the spectral weight of the elementary triplet, formally represented by the diagram in 
Fig. 6a, and 



S2{q) 



2Y,T2{k,q)^^{q,k) + 

k 



J2T,{k + ^,-k + ^,q)^|J^{q,k) 



(39) 



is the contribution of the two-particle triplet bound state, represented by the diagrams in 
Fig. 6b. Note that in calculations of the spectral weights Si{q) and 5*2 (g) the total Green's 
function, including the three particle scattering, should be used, i.e the corrected quasi- 
particle residue Z = Z^^Z^^^ and the total elementary spectrum fi^ = cu^ + Suk, should 
be substituted. Along the disordered line A = there are no quantum fiuctuations (i.e. 
Vk = 0,Uk = l,Z^'^ = 1), and the dependence of Si{q) and 5*2(5) upon 6 is given by the 
three particle scattering contribution only. In the leading order in A ^ fJ', J±, we have 
calculated the spectral weights analytically: 



Siiq) 
S2{q) 



^ ■ ^ o 2 2q^2 
-sm -11 cosig ^cos q ^coso cos -) , 

2 2^ 2 J, ^ 4J? ^ J? ^ 2^ ' 



-(-cosgsm- + - 



A 



—cos- sing)l 



(40) 
(41) 



We have also performed a self-consistent numerical evaluation of Si{q) and 5*2 (g) using 
(^,^), and the results are presented in Fig. 7a and Fig. 7b. For comparison we have also 
plotted Si{q) obtained by the dimer series expansion technique The structure factor 
of the elementary triplet is in excellent agreement with the dimer series results even for 
6 = 0.2. From Fig. 7a, which shows the case A = 0=1 — 5 — 2a, it is clear that for small a 
(and large 6), the structure factor of the bound state triplet is much smaller than the one 
for the elementary triplet. However as 6 decreases the spectral weight of the bound state 
triplet increases and eventually can even become equal to the one-particle contribution. 
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V. SUMMARY AND CONCLUSIONS. 



In summary, by using self-consistent diagrammatic analysis, we have calculated the one- 
particle and the two-particle excitation spectra of the dimerized and frustrated antiferro- 
magnetic Heisenberg chain. The quasiparticle excitations are described as a dilute, strongly- 
correlated Bose gas of triplets. The important parameter, which controls the validity of the 
diagrammatic expansion and allows us to re-sum effectively the most important diagrams 
is the density of triplets. The latter quantity is quite small for a wide range of parameters, 
making the technique, presented in this paper, effective for dimerization 6 > 0.1. For smaller 
dimerization the triplets become non-local objects and a more adequate picture would be 
one of loosely bound spinous. We emphasize however that the region of applicability of 
the presented technique is, perhaps even surprisingly, quite large. The reason for this can 
be traced to the fact that the quasiparticle weight decreases rather slowly as 6 decreases. 
Indeed, from Fig. 7a one can see that even for the quite small value 6 = 0.1 the residue is 
around 0.3 (the maximum is chosen to be 0.5). Even though this is true only in the special 
case without quantum fluctuations, it appears to be qualitatively correct even in the general 
case. Therefore the triplets are well localized objects with large quasiparticle residue even 
for rather small values of 6. We have found that our results for the spectrum are in excellent 
agreement with the numerical data, obtained by exact diagonalizations and dimer series 
expansions, for S ~ 0.2 and higher. 

We have presented a detailed diagrammatic analysis of the collective excitations and 
have found that singlet and triplet bound states of two elementary triplets generally exist 
below the two-particle continuum. In the absence of frustration there are no bound states 
in the vicinity of g = 0, while both singlet and triplet ones exist at q = it/2. Finite 
frustration increases the binding energies of both excitations, making the singlet split off 
from the continuum at g = for any a ^ 0, while a finite a is required for the triplet bound 
state to appear. We have also calculated the contribution of the elementary triplet and 
the triplet bound state into the structure factor, which is directly measurable in inelastic 
neutron scattering experiments. We have found that for small a and large 6 the bound state 
contribution is quite small, however it grows with decreasing 6 and can become equal to the 
quasiparticle contribution. 

The technique used in this work is quite similar in spirit to the dimer series expansion, 
the major difference being that while we re-sum only the most important classes of diagrams, 
the dimer series contains all of them, but only to a certain (finite) order. Let us mention 
however that at present the dimer series expansion technique has not yet been extended to 
calculate two-particle properties, such as bound state energies and their contribution into 
the structure factor ||32|. This fact makes it worthwhile, in our view, to develop and further 
refine the technique, presented in this work, especially concerning the two-particle properties 
of quantum spin models with dimerization, of which the AHC is a particular example. 
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FIGURES 



FIG. 1. Diagrams representing (second-order) corrections to the one-particle normal 
self-energy due to the three-particle scattering H^. 

FIG. 2. The elementary triplet gap, A = VI^^q/J versus a for 5 = 0.2,0.4. The squares 
represent Exact Diagonalization data |19|,20|. 



FIG. 3. Diagrams representing the three-particle scattering contributions to the two-particle 
bound states. a),b),c) for the singlet bound state, and a),b),c),d),e) for the triplet bound state. 

FIG. 4. Energy spectrum, 0,q/J, for a) 6 = 0.2, a = 0.3 (squares are Exact Diagonalization 
data mi^), h) 6 = 0.6, a = 0.0, c) <5 = 0.6, a = 0.15, d) 6 = 0.6, a = 0.2, e) 6 = 0.8, a = 0.0, 



f) (5 = 0.8, a = 0.3. The solid lines are the elementary triplet spectrum (lower curve) and the 
lower edge of the two-triplet continuum (upper curve). The dashed and dot-dashed lines represent 
respectively the triplet, S* = 1, and singlet, S = 0, bound states. 

FIG. 5. Binding energies, e^''^ / J, of the two-particle singlet, S = 0, and triplet, S = 1, bound 
states at Q = vr/2, a = versus 5, compared with the numerical Exact Diagonalization results 
(filled circles) . The solid lines represent the calculations where the three particle scattering 
is treated in simple perturbation theory, and the dashed lines represent the self-consistent inclusion 
of i?3 into the integral equation (pO). 



FIG. 6. Diagrams for the structure factors of a) the elementary triplet, and b) the two-particle 
triplet. 

FIG. 7. Structure factors Si{q) and S2{q) of the elementary and the two-particle triplet: a) 
along the disordered line 1 = 6 + 2a (A = 0), 

b) for 5 = 0.4, a = 0.2,0.4. The solid lines are the results of the calculations using Eqs. (|38| , p9| ). 
The symbols connected by dotted lines are the dimer series expansions results for Si{q). 
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